Instability development of a viscous liquid drop impacting a smooth substrate 
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We study the instability development during a viscous liquid drop impacting a smooth substrate, using high 
speed photography. The onset time of the instability highly depends on the surrounding air pressure and the 
liquid viscosity: it decreases with air pressure with the power of minus two, and increases linearly with the liquid 
viscosity. From the real-time dynamics measurements, we construct a model which compares the destabilizing 
stress from air with the stabilizing stress from liquid viscosity. Under this model, our experimental results 
indicate that at the instability onset time, the two stresses balance each other. This model also illustrates the 
different mechanisms for the inviscid and viscous regimes previously observed: the inviscid regime is stabilized 
by the surface tension and the viscous regime is stabilized by the liquid viscosity. 

PACS numbers: 47.20.-k, 47.20.Ma, 47.55.Ca, 47.55.dr 
Keywords: impact, instability, air pressure, viscosity 



The phenomenon of a liquid drop hitting a solid surface is 
ubiquitous: it occurs whenever the very first rain drop reaches 
the ground or when we spill coffee onto the floor. Liquid- solid 
impact has been extensively studied due to its broad applica- 
tions in many industrial processes, such as ink-jet printing, 
surface coating, combustion of liquid fuel, plasma spraying, 
and pesticide application! 1]. It may seem obvious that the im- 
pact outcomes should be determined by either the liquid or 
the solid properties! 2- 7], however, recent studies surprisingly 
revealed the crucial role of the surrounding atmosphere: re- 
ducing air pressure can completely suppress the liquid drop 
splashing on a smooth substrate |8l|9|], and the compressibility 
of the surrounding air is demonstrated to be important ifTol [Till . 
This unexpected discovery brings a completely new effect, 
the air effect, into the impact problem. To fully understand 
this new effect, therefore, it is essential to clarify the interac- 
tions between air and the fundamental liquid properties, such 
as surface tension and viscosity. Previous study has shown 
that the competition between the air effect and the liquid sur- 
face tension determines the impact outcomes of inviscid liq- 
uid drops! 8]. However, there has been very limited study on 
the interaction between air and the liquid viscosity, although 
the liquid viscosity itself has been broadly tested!!, B El 
and the entrapment of air bubbles in viscous drops was illus- 
trated recently lfl3l [Till . As a result, the relationship between 
surrounding air and the liquid viscosity is still missing, and 
the understanding on the liquid- solid impacts, especially the 
newly discovered air effect, remains incomplete. 

In this paper, we systematically study the interaction be- 
tween air and the liquid viscosity by varying both the sur- 
rounding air pressure and the liquid viscosity, for the impacts 
of viscous liquid drops on a smooth substrate. With high speed 
photography, we find that the instability produced by an im- 
pact highly depends on the air pressure and the liquid viscos- 
ity: the onset time of the instability decreases with air pressure 
with the power law of minus two, while it increases linearly 
with the liquid viscosity. From the real-time liquid motion 
measurements, we construct a simple model that compares the 
destabilizing stress from air with the stabilizing stress from 
the liquid viscous stress. The experimental results support the 
picture that the two stresses balance each other at the insta- 
bility onset time. This model also predicts the existence of a 



threshold viscosity, above which the system is stabilized by 
the liquid viscosity, and below which it is stabilized by the 
surface tension. This prediction quantitatively agrees with the 
previous experiment |9|]. 

We perform all the experiments inside a transparent vac- 
uum chamber whose pressure can be continuously varied from 
lkPa to 102kPa(atmospheric pressure). We also indepen- 
dently vary the liquid viscosity by using silicone oils of very 
close densities (0.92 ~ 0.94 g cm -3 ) and surface tensions 
(19.7 ~ 20.5 mN m _1 ) but different dynamic viscosities 
(4.65 ~ 13.2 mPa s). We note that all our liquids wet the sub- 
strate completely thus the wetting conditions are kept the same 
for all the impacts. To make sure that identical impact condi- 
tions are achieved each time, we release reproducible liquid 
drops of diameter d = 3.1 ± 0.1mm from a fixed height, and 
all the liquid drops impact a smooth and dry glass substrate at 
the velocity Vb = 4.03 ± 0.05 m s _1 . The impacts are subse- 
quently recorded by a high speed camera at the frame rate of 
47,000 frames per second. 

We probe the air-liquid interaction by inspecting the insta- 
bility development during the impact: under high-speed pho- 
tography, the impact produces a thin liquid film expanding 
radially along the substrate. This liquid film is stable initially, 
however, a small rim shows up around the edge at a certain 
moment, and subsequently develops into larger and larger un- 
dulations(See Fig.l left column). We believe the appearance 
of the rim indicates the transition from a stable system into 
an unstable one, and define the moment of the rim appear- 
ance as the instability onset time, t on . For example, an instant 
very close to t on is shown in the third image of Fig.l left col- 
umn. This instability onset time, t on , measures how fast the 
system goes unstable: the smaller it is, the faster the system 
becomes unstable. Interestingly, t on has a strong dependence 
on the surrounding air pressure, P. The two columns in Fig.l 
show two almost identical impacts, with only different P: At 
P = 40kPa(left column), instabilities show up in the third im- 
age; while at a higher pressure, P = 63kPa(right column), 
they appear at a much earlier time in the second image. 

By performing similar experiments under different air pres- 
sures, we systematically measure the instability onset time, 
t on , with respect to the pressure, P. We find that t on de- 
creases monotonically with P, as shown in Fig. 2. Intuitively, 
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FIG. 1. Instability development under different pressures. The liquid 
drop has diameter d — 3.1 =b 0.1mm, dynamic viscosity /i = 6.71 =b 
0.02 mPa s and impact velocity V = 4.03 ± 0.05 m s -1 . The 
left column shows an impact under the air pressure P — 40kPa. 
The impact is initially stable, but instability shows up from the third 
image. The right column shows an identical impact under higher 
pressure, P — 63 kPa. The instability appears at a much earlier time 
from the second image. 

this implies that more air leads to earlier instability appear- 
ance, thus air acts to destabilize the system, consistent with 
previous finding sJH). To test the interaction between air pres- 
sure and liquid viscosity, we perform the same t on vs. P mea- 
surements with silicone oils of very similar mass density and 
surface tension, but different dynamic viscosities, as plotted 
by the different symbols in Fig. 2. From bottom to top, the 
four curves correspond to increasing dynamic viscosities: p 
= 4.650), 6.7(o), 9.3(A), and 13.2(x) mPa s. Intriguingly, 
all the data can be excellently fitted by a simple functional 
form: t on = A/P 2 + to, with A and to the fitting parameters, 
to has typical values between 0.03 to 0.09ms, much smaller 
than most t on values. However, it is still larger than our time 
resolution(0.02ms) and can not be explained as measurement 
errors. One possibility is that the system actually becomes un- 
stable slightly earlier than the measured t on , but the instability 
features at that moment are too tiny to visualize. The pre- 
factor, A, increases with the viscosity, p, as illustrated by the 
higher locations for larger viscosity liquids. This result can be 
intuitively understood: the larger the viscosity, the more sta- 
ble the system is, and the later the instability shows up. Lim- 
ited by experimental conditions, each data set only has the 
dynamic range of about one decade in time and pressure, but 
it is nevertheless impressive that one simple functional form 
fits all the curves nicely. 

Together these data demonstrate that the instability devel- 
opment depends on both P and p, but they play opposite 
roles: P acts to destabilize the interface since higher P leads 
to faster growth of the instability; while p favors stabilizing 
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FIG. 2. (Color online) The instability onset time, t on , v.s. P for 
liquids of different viscosities. From bottom to top, the four curves 
correspond to increasing viscosities: \i = 4.65(«), 6.7(o), 9.3(A) and 
13.2(x) mPa s. All the curves can be fitted by a universal functional 
form: t on = A/P 2 + to. to ranges from 0.03ms to 0.09ms, much 
smaller than most t on values. The pre-factor A increases with p, 
as demonstrated by the higher locations of the liquids with larger p. 
Limited by the experimental condition, each data set has only about 
one decade in x and y directions. But it is nonetheless impressive 
that one simple functional form can fit all the data sets well. 



the interface as higher p slows down the instability growth. 
To quantitatively understand the effects of P and p, we in- 
spect their corresponding stresses: at the edge of the ex- 
panding liquid film, air pressure applies a destabilizing stress, 
Yjq ~ p G • Cq • V e 1 8]; and the liquid viscosity produces a 
stabilizing stress, ~ pV e /d. Here pc is the density of 
the surrounding gas, Cg is the speed of the sound in the gas, 
V e is the liquid disc expanding velocity, and d is the liquid 
film thickness measured at the edge. Cg enters the prob lem 
because previous experiments (si [2] and simulations ill (A [Tin 
suggest that the compressibility of the surrounding air is im- 
portant. 

Since V e and d vary with time, so do Eg and E M . There- 
fore a careful examination on their time dependence could 
provide valuable insight for the instability development. We 
can directly measure r(t) and d(t) from high-speed photogra- 
phy, as illustrated in Fig. 3 upper inset. V e can be obtained 
by taking the time derivative of r(t). Our measurements 
show that r(t) ex y/i, consistent with previous studies, thus 
V e = dr/dt cx 1 / y/i. This time dependence keeps valid for 
most of the expanding period, within which all our measure- 
ments are performed. Moreover, we can directly measure the 
thickness of the liquid film, d, with respect to t, as plotted in 
the main panel of Fig. 3. Because the small values of d ap- 
proach the single pixel level of our camera, the data are quite 
discrete; nonetheless they are consistent with the fit: d ~ y/vi, 
with v = p/pl being the liquid kinematic viscosity. This 
shows that d is determined by the boundary layer thickness, 
y/ut. 



From the real-time dynamics, we derive the time depen- 
dence of the stresses: The destabilizing stress, T>g ^ Pg • 
Cg -V e cx 1 / y/t, decreases with t with the power of — - ; 



3 



CM 
O 



C\J 

o 




_l I I I I 



0.1 



t(ms) 



FIG. 3. Direct measurement of the thickness d vs. time t. The impact 
is by a liquid drop of /x=4.65mPa s and Vo=4.03±0.01m/s. The inset 
shows a typical snapshot from which d is measured: d is the liquid 
film thickness measured at the edge. Main panel shows the measured 
d(t). Because d is quite small, the four discrete values correspond to 
one, two, three and four pixels of our camera. The fit is: d — 1 .9y/vi, 
indicating that d is determined by the boundary layer thickness: \fvt. 

while the stabilizing stress, E M ~ pV e /d oc 1/t, depends on 
t with the power of —1. Clearly, when t is small, E M ^> Eg, 
and the stabilizing stress dominates the destabilizing stress. 
This implies that the system should be stable initially, as we 
have observed. As t increases, however, E M decreases much 
faster than Eg and a crossover should occur at a certain time. 
After this crossover time, Eg becomes dominant and the sys- 
tem will go unstable. The experiments are consistent with this 
picture: all the impacts are indeed stable initially and become 
unstable after the instability onset time, t on . Therefore t on 
naturally corresponds to the crossover time at which the two 
stresses balance each other: 
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Plugging in the relations: pc oc P and d ex ^/Jd, with V e 
canceling each other on both sides and Cq being a constant 
independent of P, we reach the expression: 



t oc-^ 



(2) 



This expression successfully explains the two main features 
observed in Fig. 2.: (1) t on — to ex 1/P 2 and (2) the pre- 
f actor of this dependence, A, increases with /i. Moreover, 
Eq. 2 further predicts that A should increase linearly with \i. 
To test this prediction, we find A for each viscosity in Fig. 2 
from the best fit(the solid curves in Fig. 2), and plot A as the 
function of /i in Fig. 4. Indeed, a very nice linear dependence is 
observed but the line does not go through the origin; instead, it 
intercepts the x-axis at the finite viscosity value, /io = 3.4mPa 
s. 

What is the physical meaning of /io? To answer this ques- 
tion, we need to understand the impacts by the inviscid liq- 
uids with fi < jjLQ. Previous study showed that for an inviscid 
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FIG. 4. The pre-factor, A, vs. liquid viscosity, /i, for the curves 
shown in Fig. 2. The pre-factors are obtained from the best fits in Fig. 
2. A varies linearly with \i and intercepts the x-axis at /io = 3.4mPa 
S. /xo agrees with the threshold viscosity separating the inviscid and 
viscous regimes observed in previous experimentl3] . 

liquid drop impacting on a smooth surface, the destabilizing 
stress is the same as the current viscous case, Eg ~ PgCgV g 
(SO. However, the stabilizing stress, E^, is quite different. 
E l comes from the liquid surface tension, and is typically es- 
timated as the surface tension coefficient, a, divided by the 
liquid film thickness, d: E^ ~ cr/d | 8]. Therefore, we pro- 
pose that the complete stabilizing effect for an impact should 
include both the surface tension component, E^, and the vis- 
cosity component, E M . When the viscosity is small, E^ dom- 
inates E M , and we get typical inviscid behavior fl5ll . How- 
ever, when /i exceeds a certain threshold value, the viscous 
stress E M will become the major stabilizing factor, and we get 
the currently observed viscous behavior. Therefore /io natu- 
rally corresponds to this threshold viscosity which determines 
whether the inviscid or the viscous model should be used. 
We note that /io should depend on detailed impact conditions 
such as the impact velocity, surface tension and wetting con- 
ditions. Previous experiments with similar impact conditions 
already confirmed that two impact regimes exist when /i is 
varied, and the transition from the inviscid regime to the vis- 
cous regime is close to /io(see ref. |9|] Fig.5). This provides 
a strong experimental evidence for the physical meaning of 
/io. Moreover, our picture not only explains the meaning of 
/io, it also demonstrates the main difference between the two 
impact regimes: the inviscid regime is stabilized by the sur- 
face tension and the viscous regime is stabilized by the liquid 
viscosity. 

We propose that in the viscous regime, the stabilizing stress 
is mainly from the viscous stress, E M ~ pV e /d, and construct 
a model which compares E M with the destabilizing stress, 
~ PgCgV g . By assuming that Eg* and E M balance each 
other at the instability onset time, t on (Eq.l), we successfully 
explain the dependence of t on on P and /i: t on — to = Aj P 2 
and Aocfi—fiQ, with /io the threshold viscosity separating the 
inviscid and the viscous regimes. However, the most critical 
criterion, whether Eg and E M are indeed comparable at t on , 
remains to be verified. To test it quantitatively, we measure 
the ratio between the two stresses, Eg/E m ~ pcCcd/fi, at 
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FIG. 5. (Color online) The ratio between the destabilizing and the 
stabilizing stresses, Eg/E m , measured at t = t on , for various pres- 
sures and viscosities. All experiments are done under almost iden- 
tical impact conditions, with only the pressure being varied. Differ- 
ent symbols represent liquids of different viscosities: p = 4.650), 
6.7(o), 9.3(A), and 13.2(x) mPa s. The ratio, £ G /E M - p G C G d/p, 
is computed from direct measurements: p G is calculated from P, and 
d is from the best fit to the high-speed images at the time £ on . With- 
out any fitting parameter, all the ratios are within the narrow range 
between 3 and 4, confirming that Eg and E M are comparable at the 
time ton- 

the moment t on . This ratio is tested for various pressures and 
viscosities, as plotted in Fig. 5. All experiments are done un- 
der almost identical impact conditions, with only the pressure 
being varied. Different symbols represent liquids of different 
viscosities: p = 4.65(«), 6.7(o), 9.3(A), and 13.2(x) mPa s. 
For each impact, we obtain d value at t on from the high-speed 
photography measurements Blfjll . The air density pc is directly 



computed from the pressure P. The speed of sound in air at 
room temperature (20°C), Cq = 343m s _1 , is a constant in- 
dependent of P. Plugging in all the values, we obtain the ra- 
tio, Eg* /E m , as plotted in Fig. 5. Without any fitting parameter, 
most data points collapse to the narrow range between 3 and 
4. These values prove that Eg and E M are indeed comparable 
at the time t on , as our model predicts. 

We study the interaction between the air pressure and the 
liquid viscosity for the impact of a liquid drop on a smooth 
substrate. For viscous liquids, the impact is stabilized by the 
viscous stress, E M ~ I^V e /d, whose competition with the 
destabilizing stress determines when the system becomes un- 
stable. By contrast, for inviscid liquids, the stabilizing stress 
comes from the surface tension, E^ ~ a/d. Interestingly, by 
inspecting the two different stabilizing stresses, we find that 
the liquid viscosity plays opposite roles in them. For E^ in 
the inviscid regime, we have E^ ~ a/d ~ a / yfvi ex 1 / ^fp. 
Here larger p leads to larger d and smaller E^, thus more vis- 
cous liquids are less stable. However, in the viscous regime, 
we have E M ~ l^V e /d cx y/jH jfTvh . Now increasing p will 
increase E M and make the system more stable. This non- 
monotonic behavior was already observed by previous exper- 
iments(see ref. |9|] Fig. 5) and now can be fully understood. 
In summary, our study shows that the interplay between air 
and liquid viscosity is crucial in determining the outcomes of 
liquid-solid impacts. The viscosity plays different roles in dif- 
ferent regimes, and the simple intuition that a more viscous 
liquid is more stable during an impact is not always valid. 
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